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Abstract 

We isolate and study the transformation of the intrinsic spin of Dirac par¬ 
ticles as they propagate along timelike geodesics in Kerr geometry. Reference 
frames play a crucial role in the definition and measurement of the intrinsic 
spin of test particles. We show how observers located in the outer geometry 
of Kerr black holes may exploit the symmetries of the geometry to set up 
reference frames using purely geometric, locally-available information. Armed 
with these geometrically-defined reference frames, we obtain a closed-form ex¬ 
pression for the geometrically-induced spin precession of Dirac particles in the 
outer geometry of Kerr black holes. We show that the spin of Dirac parti¬ 
cles does not precess on the equatorial place of Kerr geometry; and hence, in 
Schwarschild geometry. 
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1 Introduction 

Our goal in this paper is to isolate and study the transformation of the intrinsic spin 
of Dirac particles as they propagate along timelike geodesics in Kerr geometry. The 
motivation comes from the problem of observers located in the outer geometry of 
Kerr black holes trying to communicate quantum information by exchanging polar¬ 
ized particles. It was shown in [FKP14] how observers located in the outer geometry 
of Kerr black holes may exploit the symmetries of the geometry in order to com¬ 
municate information by exchanging polarized photons. In the present paper, we 
consider the problem of observers in Kerr geometry trying to communicate quantum 
information by exchanging massive spin-| particles instead. In the massive spin-| 
case, one encounters two problems that are not present in the photonic case covered 
in |FKP14j . which we now describe. 

Suppose Alice encodes information in the spin of a massive spin-f particle and 
sends it to Bob. First, she needs to ensure that the particle’s trajectory will intersect 
Bob’s worldline. Spin-f particles are described by spinor fields that obey the Dirac 
equation. It is not immediate how to relate a given solution to the Dirac equation, a 
spinor field, to a timelike geodesic along which the test particle propagates. Second, 
whereas a photon’s polarization 4-vector is parallel propagated and remains orthog¬ 
onal to its 4-velocity, there is no such propagation law for the spin vector of a Dirac 
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particle. Following |Aud81j . we address these two problems by using the semiclas- 
sical ansatz for the Dirac equation, which allows us to recover a timelike geodesic 
along which the test particle propagates and derive a propagation law for the spin 
vector. 

Terashima and Ueda’s seminal paper Tb() 1! outlined a strategy for evaluating 
the spin precession induced by the motion of spin-| particles in a curved space- 
time; a strategy that has since been followed by numerous authors [AJK09 . AF12L 
IPTW121ILanl2irSA101lRPGCllj . The strategy outlined in |TU04] is to calculate the 
Wigner rotation induced by the instantaneous local Lorentz transformation relating 
the particle’s 4-momentum at nearby events along the particle’s worldline. For a 
particle moving along a timelike geodesic, the strategy assumes that the precession 
of the intrinsic spin of the particle is determined solely by the rotation of the ob¬ 
server’s frame with respect to a frame that is parallel propagated along the particle’s 
geodesic worldline0 jRPGClll ISAlOj follow the strategy of [TU04 J to obtain the 
spin precession of spin-| particles on circular orbits confined to the equatorial plane 
of Kerr-Newman geometry. |Lanl2| also follows the strategy of [ TU04j to obtain the 
spin precession of spin-^ particles on circular and radially-infalling geodesic orbits 
confined to the equatorial plane of Kerr geometry. 

We take a very different approach from |TU04j . We work directly with the Dirac 
equation and, following |Aud81] , use the decomposition of the Gordon decomposition 
of the Dirac current to first define the spin vector of a Dirac particle. Then we intro¬ 
duce the semiclassical ansatz for the Dirac equation and thereby recover the geodesic 
along which the Dirac particle propagates; a strategy first proposed in [Pau32j. We 
extend the results of [Aud81] and show that the spin vector is parallel propagated 
along the aforementioned geodesic to 0(h 1 2 ) in the semiclassical ansatz (Theorem 
15.2p . We develop a new method of constructing a reference frame on purely geomet¬ 
ric criteria that allows us to obtain the proper time-dependent rotation of the spin 
vector in a coordinate independent manner. Our expression for the geometrically- 
induced precession of the spin vector is valid for a Dirac particle propagating along an 
arbitrary timelike geodesic in the outer geometry of a Kerr black hole. We also obtain 
an expression for the spherical curvature of the curve traced out by the spin vector, 
which allows us to analyze the dynamical behaivour of the spin vector with a single 
invariant function. Even though our approach is quite different from [TIJ04] . our 
qualitative result that the spin vector is parallel propagated along a Dirac particle’s 
geodesic worldline agrees with the assumption underlying their strategy. 

The rest of this paper is organized as follows. In Section [2] we recall some of 

1 If the particle is accelerated, there is an additional term that arises from boosting the 4- 

momentum along the worldline of the particle. 
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the salient geometric properties of the Kerr metric that will be used in this paper; 
including the equations of motion for timelike geodesics and parallel propagated 
frames along timelike geodesics. In Section [3J we introduce the Dirac equation, 
recall the Gordon decomposition of the Dirac current and define the spin vector of 
a Dirac particle. Next, in Section [4] we introduce the semiclassical ansatz for the 
Dirac equation and thereby recover the classical trajectory of a Dirac particle. We 
prove a propagation law for the spin vector in Kerr geometry in Section [5] Then, we 
construct a reference frame on purely geometric criteria and define the geometrically- 
induced precession of the spin vector in Section [ 6 j In Section [Tj we provide an explicit 
expression for the geometrically-induced precession of the spin vector and obtain a 
spherical curvature invariant for the curve traced out by the spin vector. We conclude 
with a discussion of the findings in Section [HI 

Remark 1.1 (Notation). We shall reserve lower case Latin indices, a, b, c,..., for ar¬ 
bitrary orthonormal frames in which the metric takes the form g ab := diag(l, - 1 ,- 1 , 
— 1); bracketed Latin indices, (a), ( 6 ), (c),..., for orthonormal frames that are par¬ 
allel propagated along a timelike geodesic and in which the metric takes the form 
:= diag(l, —1, —1, —1); hatted Latin indices, i, j, k ,..., for the spacelike com¬ 
ponents in an orthonormal frame; unhatted Latin indices i,j , k,..., and Greek in¬ 
dices, a, (3, 7 ,..., for spacetime coordinates in which the line element takes the form 
ds 2 = g a pdx a dxP. We will, on occasion, use index-free notation as follows. Given 
a metric ( g y) and vector field X = the 1-form dual to X will be denoted 

by X b , whose components are given by (X b )i := Similarly, given a 1-form 

cj := u>idx\ the vector held dual to u will be denoted by Cod, whose components 
are given by (u^)* := c cjg l K We shall sometimes find it convenient to use semicolons 
to denote covariant derivatives, e.g., T ;Q , := V a T, whereas commas will denote or¬ 
dinary partial derivatives, e.g. f, a „/• We shall work throughout in natural 
units: G = c = 1. 

2 Kerr geometry 

Remark 2.1. We collect in this section some well-known facts about the Kerr metric. 
Further details can be found in |FKP14j . 

In Boyer-Lindquist coordinates ( x l ) = (t,r, with —oo <t< +oo, r + < r < 
Too, 0 <r 9 < 7 r, 0<(^< 27r, the Kerr metric takes the form 

ds 2 = — (dt — asm 2 ddip) 2 — ~rdr 2 — 'Edd 2 — v (adt — (r 2 T a 2 ) dp) 2 , (2.1) 

2-j ZA 2.J 
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with 


£(r, P) := r 2 + a 2 cos 2 P, A(r) := r 2 — 2Mr + a 2 . (2.2) 

The parameters M > 0 and a > 0 correspond respectively to the mass and angular 
momentum per unit mass of the black hole, as measured from infinity. We shall only 
be considering the non-extreme case M > a > 0, which implies that the function 
A(r) has two distinct zeros, 

r± = M± VM 2 - a 2 . (2.3) 

Moreover, we shall restrict our attention to the region r > r + , which describes the 
geometry outside the event horizon of the black hole. 

The Kerr metric admits a two-parameter Abelian isometry group generated by 
the pair of commuting Killing vector fields d t and d^. The Kerr metric also admits 
a discrete subgroup isomorphic to Z 2 generated by the involutive isometry 


( t , r, P, <p) ^ {-t, r, P, -ip). (2.4) 

We shall denote by L the differential of the isometry (12.41) . The Weyl conformal 
curvature tensor of the Kerr solution is of Petrov type D, meaning that it admits a 
pair of repeated principal null directions, each of which is defined up to multiplication 
by a non-zero scalar function. We eliminate the scaling freedom we would have 
otherwise had in defining a null coframe adapted to the principal null directions of 
the Weyl tensor as follows. 

Definition 2.2 (Symmetric frame). Our null coframe is chosen such that 

LP 1 = -P 2 , L-d 2 = -P 1 , LP 3 = -tf 4 , Ltf 4 = -tf 3 . (2.5) 

We refer to this frame as the symmetric null coframe. It is given in Boyer-Lindquist 
coordinates by 

■d 2 


a/2SA 

1 

V2SA 

1 


(A dt + Sdr — a sin 2 PA dp) , 
(A dt — Yidr — a sin 2 PA dp) , 


V2S 

1 

V2S 


((r 2 + a 2 ) sin i9d<p — iSdP — a sin Pdf) , 
(P + O sin tfdip + iT.d'd — a sin ddt j . 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 

(2.9) 
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The orthonormal symmetric coframe (u; 0 , uj 1 , lo 2 , a; 3 ) corresponding to the sym¬ 
metric null coframe (i? 1 , i? 2 , i? 3 , f? 4 ) is then defined by 


+ tf 2 ), w 1 = - O 1 ) V = -^(O 3 + <3 4 ), u. 3 = f^(0 3 - 

( 2 . 10 ) 

The principal null directions of the Weyl tensor, with scale factors chosen according 
to the requirement (12. 5 j) . are given by 


£ = eS- = 


and 


dx i V2SA 
d 1 


n = n 


ox d f-r- d d \ 

+ a)- + N / A- + 0 —), 


+ “ 2) Ft ~ + %)' 


( 2 . 11 ) 


( 2 . 12 ) 


dx i a/2SA 

We now define observers in terms of the principal null directions of the Weyl tensor 
and the involution L. 


Definition 2.3 (Carter observers). The vector field, 

°'- = 75 {e + n)= vis ( (r2 + " 2) s + %) ’ (213) 

where £ and n are given by (12. Ill) and (I2.12jh is timelike and future-pointing, and 
identifies a family of observers, that we call Carter observers, whose 4-velocities are 
symmetric linear combinations of the principal null directions of the Weyl tensor. 

Carter observers exist everywhere outside the event horizon including the region 
between the event horizon and the stationary limit surface where the stationary 
Killing field £ = d t becomes null. Their angular velocity is r2 f a 2 , which is exactly 
the angular velocity of the event horizon; both as measured at infinity. Therefore, 
this class of observers is uniquely suited to analyze the behaviour of test particles 
near the horizon. 

We choose the observers’ frames to be dual to the symmetric coframe defined by 
(I2.10K . These frames can constructed using locally-available geometric data. Specif¬ 
ically, the construction of the symmetric frame only requires knowledge of the prin¬ 
cipal null directions of the Weyl tensor. 

In addition to its two-parameter Abelian group of isometries and the involutive 
isometry, the Kerr metric posesses further symmetries that are “hidden” in the sense 
that they cannot be represented by Killing vector fields. The existence of these hidden 
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symmetries is closely tied to the fact that all the known massless and massive wave 
equations are separable. The geometric object that generates all these additional 
symmetries is a rank two Killing-Yano tensor , that is, a (0, 2) skew-symmetric tensor 
( fij) satisfying the Killing-Yano equation, 

^(ifj)k = 0. (2.14) 

In Boyer-Lindquist coordinates and in the symmetric orthonormal coframe, 

f = — acos'duj 0 A us 1 + ruj 2 A u> 3 , (2.15) 


is a Killing-Yano 2-form. The symmetric (0,2)-tensor defined by 


K,j 



(2.16) 


satisfies the Killing equation 


^ ( i^jk ) 0 ? 


(2.17) 


and therefore gives rise to a quadratic first integral for the geodesic flow in Kerr 
geometry first discovered by [ Car68] . 


k = KijlTW, 


(2.18) 


where (U l ) is the 4-velocity. The quadratic first integral defined by equation (12.1 8ft 
exists in addition to the two linear first integrals arising from the presence of the two 
commuting Killing vector fields d t and and therefore reduces the integration of 
the geodesic flow to quadratures. 

The equations of motion for timelike geodesics are given in Boyer-Lindquist co¬ 
ordinates by 



£A i 

= ((r 2 + a 2 ) 2 - a 2 A sin 2 E - 2Mra<h, 

(2.19) 



_ r _. A — a 2 sin 2 ft T 

= 2 MraE -|-^—-$, 

sin" v 

(2.20) 


S 2 r 2 

= R(r), 

(2.21) 


S 2 ^ 2 

= ew, 

(2.22) 

where 


R(r) := P 2 — A (ft + r 2 ) , 

(2.23) 



0 ($) := k — acos 2, d — D 2 , 

(2.24) 


7 







and 


(2.25) 

(2.26) 


P(r) := E (r 2 + a 2 ) — a$, 

$ 

D(t?) := aE sill'd- 

sinw 

In equations (12. 19112.22]) . E := p t is the conserved energy, $ := — p v is the conserved 
angular momentum, and the affine parameter is chosen such that the 4-velocity ( U l ) 
has unit norm, gijU l U J = +1. The 4-velocity (U l ) of an arbitrary timelike geodesic 
is given in the symmetric coframe by 


U := 


1 


p VR 

v/A’v/A’ 



(2.27) 


In order to obtain an orthonormal frame that is parallel propagated along U, 
we recall Marck’s elegant construction |Mar83j . 

Note first that U is parallel propagated along itself so that we can set 


V : = u °■ < 2 - 28 ) 

The Killing-Yano 2-form /, defined by equation (12. 14[) and given in the symmetric 
frame by equation (j2. 15(1 . gives rise to a spacelike vector ( L a ), 

L a ■= f\U b , (2.29) 


that is parallel propagated along U as a direct consequence of the Killing-Yano 
equation. We normalize L , defined by equation (I2.29[) . to obtain a unit-norm vector 
given by 



a cos fly/R acos'dP 


v/A 


VA 


, ry/Q, —rD 


(2.30) 


In order to obtain the remaining two vectors of the parallel propagated frame, we 
start with a basis 


L (i) 

L (2) 


1 / wryjR zurF acos'd-y/© acos'dD 

V \/A ’\/A’ w ’ w 
1 / wF vjy/R D i/©\ 

\/E V \/A ’ v'A w 5 w J ’ 


(2.31) 

(2.32) 


where 


k — a 2 cos 2 'd 


w : = 


r 2 + k 


(2.33) 





























Then we solve the ODE that governs a proper time dependent rotation angle <f>(r) 
such that the two vectors 


L{ i) = cos <E>(t ) A(i) — sin <3>(r) A (2 ), (2.34) 

L( 2 ) = sin$(r)A ( i) + cos$(t)A( 2 ), (2.35) 


are parallel propagated along U. The solution, originally obtained in |Mar83] . is 
given by 


dr 


1 

K2 


P 


r 2 + k 


a sin id 


k — a 2 cos 2 id 


(2.36) 


where r denotes proper time. The proper time dependent angle $ may be obtained 
by separation of variables, 


4>(r, id) := F(r) + G(-&), 


(2.37) 


where 


F(r) 


:= K 2 


P dr 
r 2 + K y/R ’ 


G(ld) = CLK 2 


sin i93 did 
k — a 2 cos 2 id i/O 


(2.38) 


Theorem 2.4 (Marck). Let 7 be a timelike geodesic parametrized by proper time 
t, with tangent vector U := ^ 7 . Then, the orthonormal frame, L( a ), defined by 
equations H2.28\) . H2.34\)> A 2. 3,51) . and H2.30\) . is parallel propagated along 7 . 


3 The Dirac equation 

The Dirac equation is given by 


ih^ a W = m\ I/, (3.1) 

where T is a 4-component spinor field, m > 0 is the rest mass of the Dirac particle, 
V a is the covariant derivative, h is Planck’s constant, and 7 “ is a representation of 
the Clifford algebra Cf?i j 3 (M), 

7 a 7 6 + 7 V = 27 a T (3.2) 

where 1 is the unit element of the algebra. 

As a direct consequence of the Clifford algebra relations (13.2[) . the commutators 
of the gamma matrices, 

(3.3) 
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constitute a representation of the Lie algebra s 1(2, C) defined by the relations 


[< 7 ab , a cd ] = i (jj 


adg.bc 


r] ac a bd + rf c o ad 


v bd a ac 


)• 


(3.4) 


The generators of the Lie algebra fl3.3jb satisfying the commutation relations (13.4ft . 
will play a key role in our definition of the spin vector for Dirac particles. We shall 
be using the standard representation of the Clifford algebra given by 


7° = 


/ 0 
0 -I 


7 


0 cP \ 
-cP 0 ) 


(* = 1,2,3), 


(3.5) 


where I is the 2x2 identity matrix. The 2x2 Pauli spin matrices cP are explicitly 
by 


cr 


i . 


0 1 
10/’ 


<x 2 := 


0 i \ 

i or 


a 3 := 


1 0 
0 -1 


(3.6) 


Given a spinor field T, the adjoint spinor is defined by T := fCy 0 , where the dag¬ 
ger denotes complex conjugation and transposition. If T satisfies the Dirac equation 
(13.11) . then T satisfies the adjoint Dirac equation , 


i/iV a f , 7 a = -mf. 


(3.7) 


For a spinor field T satisfying the Dirac equation (13.ip . the quantity 


f := T 7 a T, 


(3.8) 


defines a vector field called the Dirac current. The Dirac current is conserved, that 
is, 

V a f = 0, (3.9) 

as can be seen by multiplying equation (13.1ft on the left by T, multiplying equation 
(13.71) on the right by T, and adding the two terms together to eliminate the mass 
term. 

We now turn to Gordon’s result }Gor28j . which demonstrates that the Dirac 
current decomposes into two parts which are separately conserved. 

Theorem 3.1 (Gordon decompositon of the Dirac current). The Dirac current de¬ 
composes into convection and polarization 4-currents which are separately conserved. 
More precisely, 

=--j a =j a pola r+j a con, ( 3 - 10 ) 

where 

*<•" := i (^ ■ I L- (3 - u) 


10 














is the polarization current satisfying V a jp 0 i ar = 0 and 

fccn ■■= ff,Ti - ® $; “) > < 3 ' 12 ) 

is the convection 4-current satisfying 'V a jf on = 0. 

Proof. The decomposition follows from noting that 

j = 


where we have used the Dirac equation and its conjugate repeatedly, along with the 
symmetrization and antisymmetrization of the Dirac gamma matrices. The conser¬ 
vation of the convection current follows from recalling Lichnerowicz’ identity [Lic64j . 

A = -V a V a + ^R, (3.16) 

where R is the scalar curvature and A is the Laplace operator whose action on the 
Dirac spinor and its adjoint is given by 


(T 7 VT) ;fe -(T 7 VT) ;fe (3.13) 

9 777 , _ _ 

— j a + T. b7 a 7 6 T - T 7 6 7 a ^& (3-14) 

in 

4.777 _ _ 

— f + 2 (T ; “T - TT ;a ) , (3.15) 


/ TYl \ ^ 

AT = 7 ° 7 & V a VftT = (-J T, (3.17) 

and 

AT = V a V fe T 7 V = (-J T, (3.18) 

respectively. It follows that 

(T ;a T _ TT ;Q ). a = (V“V n T) T — T (V a V a T) (3.19) 

= ii?TT - (AT) T - ii?TT + T (AT), (3.20) 

= (^-rj (-TT + TT) (3.21) 

= 0. (3.22) 


Since both the Dirac current and the convection current are conserved, the polariza¬ 
tion current, defined by equation (13.lip , must be conserved as well. □ 
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For a given Dirac spinor 'I', the convection current j“ on defines a congruence of 
timelike curves with unit tangent vector (K a ) defined by: 


K a := 


rt dt 

J con 


\ZVabjcon 


jb 

J con 


(3.23) 


We are now in a position to define the spin vector. The following definition was 
originally proposed in l AndSl . 


Definition 3.2 (Spin vector). We define the spin vector associated to the spinor T 
by 



abed 


K b 




(3.24) 


where e abcd is the anti-symmetric Levi-Civita symbol given by 


^.abed . 


+1 if abed is an even permutation of 0123, 
— 1 if abed is an odd permutation of 0123, 
0 for repeated indices. 


(3.25) 


Since the cr-matrices, defined by equation (13. 3 ll . are the generators of the Lie al¬ 
gebra (13.41) . we may interpret 'Lc 7 a6 'F as spin density. Kirsch et al. have shown that 
the spin operator defined uniquely via the Gordon decomposition corresponds to the 
Foldy-Wouthuysen mean-spin operator [KRH01] , Meanwhile, Bauke et al. have re¬ 
cently shown that, in the Lagrangian formulation, the Foldy-Wouthuysen mean-spin 
operator is the only known relativistic spin operator that commutes with the free 
Dirac Hamiltonian, has the eigenvalues ±|, and obeys the angular momentum al¬ 
gebra |BAKG14) . Similar results were derived in (CRW13) . who further show that 
the Foldy-Wouthuysen mean-spin operator is the only spin operator proposed so 
far, that has the right non-relativistic limit and does not convert positive (negative) 
energy states into negative (positive) energy states (their “charge symmetry condi¬ 
tion”). These results give us good confidence that the definition (I3.2j) . first proposed 
in |Aud8lj . of the spin vector in terms of the Gordon decomposition of the Dirac 
current is a good one. 

We now turn to the semiclassical ansatz which will allow us to recover the classical 
trajectory of a Dirac particle and derive a propagation law for the spin vector defined 
by 0323) . 
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4 Semiclassical ansatz 


The semiclassical ansatz for the Dirac equation is a formal power series expansion in 
Planck’s constant h, and given explicitly by 

OO 

T(x) = exp(iS(x)/h) (— ih) n a n (x). (4.1) 

n =0 

where S{x) is a scalar field and a n (n > 0 ) is a countable sequence of 4-component 
spinor fields. Plugging the ansatz (14.ip into the Dirac equation (13.11) and setting to 
zero the coefficients of the different powers of h, we obtain 

(7 a S )a + mI)a 0 = 0, (4.2) 

(l a S, a + ml) a n+ i = - 7 a a n - a (Vn 6 N), (4.3) 

The existence of a solution to the homogeneous equation (14. 2 j) requires 

det( 7 a S' iQ + ml) = 0. (4.4) 


Equation (14. 4 j) in equivalent to the Hamilton-Jacobi equation for timelike geodesics, 


S’ a ,S, Q = m 2 . 

(4.5) 

We define 


Pa ■ ‘S',aj 

(4.6) 

and normalize to obtain a unit-norm vector field 


U a := —p a = -—S’ a . 

(4.7) 


m m 


As a result of the Hamilton-Jacobi equation (14.51) . the integral curves of U , defined 
by (14.71) . are guaranteed to be timelike geodesics. And since equation (14.21) is the 
classical limit (fl —> 0) for the Dirac equation with the semiclassical ansatz, we 
interpret p a , given by equation (14.6p . as the 4-momentum, and f/ a , given by (14.7[) . as 
the 4-velocity of the Dirac particle described by T. 

The use of the semiclassical ansatz for the Dirac equation (14. lj) and the obser¬ 
vation that the homogeneous equation (14.21) implies the Hamilton-Jacobi equation 
for spinless particles goes back to Pauli [Pau32j . Bargmann et al. derived equa¬ 
tions describing the classical trajectories of spin-^ particles in uniform and constant 
electric and magnetic fields [BMT59) . Rubinow and Keller showed that the classi¬ 
cal equations of Bargmann et al. could be obtained from the Dirac equation using 
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the semiclassical ansatz [ RK63j . Rafenelli and Schiller obtained essentially the same 
result soon after, using the so-called “squared” Dirac equation along with the semi- 
classical ansatz |RS64j . 

We follow the more recent work by |Aud81j who showed that the spin vector, 
defined by equation (I3.24p . is parallel propagated along the particle’s trajectory to 
zeroth order in the asymptotic expansion (14.11) . We will extend the result of [Aud8lj 
and show that, in Kerr geometry, the spin vector is parallel propagated along U a to 
first-order in h. We begin our analysis with recalling some results from |Aud8l] , 

The matrix acting on a 0 in equation (j4.2[) is of rank 2, as is manifest by con¬ 
sidering a frame that is parallel propagated along the congruence U a in which the 
4-momentum takes the form p^ = (m, 0, 0, 0). The matrix acting on a 0 in equation 
(14.21) is then given by 


( 0 

0 

0 

0 

\ 

0 

0 

0 

0 


0 

0 

2m 

0 


^0 

0 

0 

2m 

J 


The general solution of the homogeneous equation (|4.2[) therefore takes the form 


0-0 — Al&Ol + /?2&02, 


(4.9) 


where the basis 4-spinors b m and b 02 are two linearly independent solutions 


An — 


E + m 
2m 


1 


0 

0 

E + m 

1 

p 5 

E+m 

' 602 “V 2m 

p X —ip * 

E+m 

p 1 +ip 2 


P 3 

. E+m . 


E+m _ 


(4.10) 


where E := p° is the energy of the Dirac particle and p a is defined by equation (14.6|) . 
The basis spinors (I4.10|) are parallel propagated along U, 

boi; a U a = 0, b 02 . a U a = 0. (4-11) 


Choosing a frame that is parallel propagated along the congruence U a in which the 
4-momentum takes the form p^ = ( m , 0, 0, 0), the spinors b 0 1 and b 02 reduce to 


' 1 ' 


1 

O 

_ 1 

0 


1 

0 

) fy )2 — 

O 

1 

O 

_1 


- 1 

O 


(4.12) 
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With the inner product on spinors defined by := xV’ = i^Xi f° r arbitrary 

4-spinors y, -0, an orthonormal basis for the 2-plane orthogonal to the fundamental 
solutions (jUmp of the homogeneous system (14.21) is given by the spinors 


611 = 


E + m 
2m 


r 1 


1.9“ 

p L —ip* 

E+m 


E+m 

pWip 2 

E + m 

P 3 

E+m 

, O12 — \ 

E-\-m 

1 

V 2m 

0 

0 


1 


(4.13) 


The spinor fields (b 0 i, b 02 , bn , 612 ) constitute a basis for 4-spinors. In a parallel prop¬ 
agated frame along the congruence U a in which the 4-momentum takes the form 
p( a ) _ ; q, 0 ), 611 and 612 take the form 


-1 

O 

_1 


1 

O 

_ 1 

0 


0 

1 

, b \ 2 — 

0 

1 

O 

_1 


1 


These spinors are parallel propagated along U as well, 


(4.14) 


bn„U“ = 0 , i>i 2 i „C/ a = 0 . 


The general solution to (14.31) for n — 1 can be written as 


di — Viboi + v 2 bo2 + Ai 6 n + X 2 b\ 2 . 


(4.15) 


(4.16) 


We are now in a position to state the result from [ AudSlj that would allow us to 
derive a propagation equation for the spin vector in Kerr geometry. 

Lemma 4.1 (Audretsch). The propagation equation for the scalar coefficients (/3 = 
j3\, j3 2 , Vi, v 2 ) of the fundamental spinors An and 602 in the general solution j4.16\ ) to 
the first-order equation o with n = 1 is given by 

A aU‘ = (4.17) 

In order to prove that the spin vector defined by (13.241) is parallel propagated to 
first-order in the asymptotic expansion (14.lj) . we will now to determine the solutions 
of the propagation equation for the scalar coefficients in Kerr geometry. 


15 





























Theorem 4.2 (Scalar coefficients in Kerr geometry). Let U be tangent to a timelike 
geodesic in Kerr geometry. The general solution to the propagation equation for the 
scalar coefficients \f.l r i\) is given by 


P(r,ti) ■ = 


a/.R( r)0($) sin$’ 


(4.18) 


where c is a constant of integration. 

Proof. We may rewrite the propagation equation (14.171) in terms of proper time r as 
follows. 


d 

dr 


In P{t) 


= 

(4.19) 

1 (R! ©' • A 

(4.20) 

= — — r + —d + cot dd , 

2 \ R 0 J 

= — ^^-ln(i?0sin^), 

2 dr 

(4.21) 


where the dot denotes prime denotes ordinary partial derivatives, and where the 
functions R(r) and ©($) are defined by equations (12.231) and (12.241) respectively. It 
follows from equation (14.2ip that the general solution to (14.17j) is given by the scalar 
field (14481) . □ 


Let 


We can easily verify that 


Pi 

P 2 


Vl 


V2 


_Cl_ 

a/ R(r)Q(i!)) sill'd’ 
_ C2 _ 

a/-R( r)©($) sill'd’ 
d\ 

a/ R(r)Q(f) sin^’ 

c?2 

a jR{r)Q{f) sin d 


f ao«o — y/PlPi + P 2 P 2 , 


(4.22) 

(4.23) 

(4.24) 

(4.25) 

(4.26) 


satisfies the propagation equation 


laU a 



(4.27) 
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Define a spinor bo such that 


ao = fb Q . 


Then, b 0 has unit norm, 


bobo — 1, 


and is parallel propagated along U, 


b 0 ;aU a = 0. 


(4.28) 

(4.29) 

(4.30) 


We now recall our definition of the spin vector given by equation (13.24ji , 


W a l £ abcdj^cd^_ 

2 TT ’ 


(4.31) 


where e abcd is the Levi-Civita symbol. In order to derive a propagation equation for 
the spin vector to first-order in h, we note that 


1 


— ~C9 ( + ~P? («0«1 — OlOo) ) + 0{K 2 ). 


vpvl; p \ p 
Therefore, the spin vector is given by 

W a = W$ + h\V° + 0(h 2 ), 


where 


and 


Wo a = - e abcd U b b 0 E cd b 0 , 


ra _abed 


W,° =-£° 
1 2 


U b 


aoai — aiaor , 1 _ \ 

- -p - o 0 a cd o 0 - -p ( a 0 a cd ai - aia cd a 0 ) 


— 7^ (bo;bbo ~ bobo-b) bocr cd bo 


The spin vector defined by (13.241) and given by (14.331) satisfies, 

W a W a = -1 + 0(h 2 ). 


(4.32) 

(4.33) 

(4.34) 


(4.35) 


(4.36) 


In order to prove that the spin vector defined by (13.241) is parallel propagated to 
first-order in the asymptotic expansion (14.ip in Kerr geometry we will show that, 
to first-order in h, the spin vector defined by (j3.24[) has constant components in the 
parallel propagated frame. 
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5 Propagation law for the spin vector 

In order to prove that the spin vector defined by (13.241) is parallel propagated to 
first-order in the formal asymptotic expansion (14. Tj) in Kerr geometry, we will need 
the following lemma. 


Lemma 5.1. For the general solution to the propagation equation for the scalar 
coefficients in Kerr geometry /3(r,D), given by we have 


A(0) = 


A(i) = - 


A(2) - - 


£_ 

2E 

A 

2E 


wR 1 ^cos 2 — gin + 


cot fl + 0'/0 


w 


(5.1) 

cos $a cos "dDfi ; - 2 — sin $0 
(5 ' 2) 

wR' ^sin$rPfi ; _ 2 + cos<I>j -|- t ^ + 0 /© ^ s in <ha cos^O/t - ^ + cos*!)©^ 

(5.3) 

7 


A(3) - 


2^2 £ 


„ R! 

( n 0 'Y 

—a costdP-h rD 

cot d + — ) 

R 

l 0 ;J 


(5.4) 


where the prime denotes ordinary derivatives. 


Proof. Let ( X , Y) := X a Y a denote the inner product for arbitrary vector fields X and 
Y. Observe that 7 (a) = (d/3,L^ a f). The proof then follows from direct computation 
using the parallel propagated frame L( a ) constructed in Section ??. 


(d/3,L (0) ) = (d/3,U) = - 9 -/3. 


(d/3, L(i)> = 


7 


2^2 £ 


m R a cos'd, 

r©P-1-1 

R w 


a 07 

cot d + — 


(d(3,L {2) ) = - 


2E 


0 


wRf H-( cot D + 


w 


&_ 


(d/3, L(3)) — 


7 


2k2S 


„ K 

( „ &y 

—a cos'dP-b rD 

cot fl + — ) 

R 



(5.5) 

(5.6) 

(5.7) 

(5.8) 

□ 


Theorem 5.2. In Kerr geometry, to first-order in h, the spin vector defined by \3.2f\) 
is parallel propagated along U. 
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Proof. Theorem 14.21 guarantees that the scalar coefficients of the spinor fields take 
the form fl4.22jl - fl4.25h for some constants of integration Ci,Ci,d\,d 2 . Plugging the 
scalar coefficients fl4.22jl - fl4.25jl into the expressions previously established for the 
spin vector, (j4.34jl and f!4.35f) . we obtain in the parallel propagated frame 


and 

where 


W n = 


Ci Ci + C 2 C 2 


0 

Ci C2 + C%C\ 

i(c*c 2 - C 2 C 1 ) 

C* 2 C 2 - Ci Ci 


Wi = P + Q + R, 


o_ if. - + 

P “tv (®o a i — affio) Wo — 1 ----—-- Wo, 

J C-^Ci + C 2 C2 


Q = 2^ (a)(0)u { a i°Tj a o - Ooc-Oi) , 


Cl Cl + C 2 C2 


0 

{d\c 2 + d 2 ci) — (c\d 2 + c 2 d\) 
i {{c 2 d\ + d*C2) — {c\d 2 + d 2 cf)) 
(c*di + d 2 c 2 ) — {d\c\ + c* 2 d 2 ) 


and 


and where 5J is given by 


R = -eWWMMi 


- * 0 * 01 ( 0 )) 


(5.9) 

(5.10) 

(5.11) 

(5.12) 

(5.13) 

(5.14) 

(5.15) 


Remark 5.3. Note that P and Q have constant components in a parallel propagated 
frame. We will now show that R , given by (15. 14j) . vanishes identically. 

We begin by evaluating (J5T 


fr();(a)fro — bobo-(a) 


{faL(a) {fa) + fafaa) (&*)) - {fafaa) {fa) + fafaa) {fa)) 

n 


2 Recall that do = fb 0 = ffifooi + /3 2 bo2 and / = + P 2 P 2 in accordance with (14.261) . (14.281) . 

and (TOR) . 
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(5.17) 


where L( a )(/3) := (d/3,L( a )). Since L( 0 )(/3) = (d/3,U) = — f/3, we have 

A;(0)A — AA;(0) = 0. 

And for the spacelike basis vectors of the parallel propagated frame L^, 

PlP* ,(t) + ^2^2,(i) ~~ PlPl,(i) _ /3|/?2,(i) 


k);(i)A> — «0°0;(i) 


bnbn-(i\ — 


/ 2 


(5.18) 


Since 


&oA<b(o) fc o = 0 and 6 0; (o) = 0, (5.19) 

at least one of the three terms, eWWWW, <5J( a ), and b 0 cr^)(c)bo, vanishes for any 
choice of indices (a), (5), (c), one of which has to be (0). Thus, R (l) must vanish for 
i = 1,2,3. 

The timclike component of R is given by 


R {0) = 


1 


2 mif 4 


(Pi 02 + @2Pi) {PudPI — Pt(l)Pl + A,(1)A> — A>,(1)A) 


+ i (P 1 P 2 — P 2 P 1 ) (Pl,( 2 )Pl — Pi,{ 2 )Pi + A,(2) A — A,( 2 ) A) 
+ (P 2 P 2 — Pi Pi) (A,(3) A ~ Pi, (3) Pi + A,(3) A ~ A,(3)A) 


(5.20) 


Note that we have a common factor between the derivatives of the scalar fields (15.2ft - 
(15.41) . Define 

z\ Alb _ Alb _ ^Ib _ ^(i) r 5 

CU ' A P 2 PI P* 2 ’ 

Therefore we may replace each term, A,(b> by C (i)Pk (for k = 1, 2, 3) in the expression 
for R^ given by (15.201) as follows. 


R w = 


C(l) (AA + AA) (AA ~ AA + AA - AA) 


2 mif 4 

+ i(( 2 ) (AA - AA) (AA* - pip ! + AA - AA) 
+ C(3) (AA - PI Pi) (AA - PI Pi + AA - AA) 
= 0. 


(5.22) 


This completes the proof of our claim that R vanishes identically. 

Since (ci,C 2 , A, A) are constant along the integral curve of U, we see that, to 
first order in h, the spin vector has constant components in the parallel propagated 
frame, and is thus parallel propagated along U. □ 
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For the sake of completeness, we note the expression of the spin vector for a Dirac 
particle in Kerr geometry to first-order in h. 

w a = W 0 a + hW? + 0(h 2 ), (5.23) 


where 


Wa = 


cjci + c * 2 c 2 


0 

c\c 2 + c 2 c\ 
i{d[c 2 - c* 2 d) 

C2C2 - cld 


and 


Wi = - 

C*Ci + c*c 2 
where c 0 := (c*di + c* 2 d 2 ) — {d\c\ + d 2 c 2 ) 


( 

0 

\ 

Co Wo + 

( d\c 2 + d 2 c\) — ( c[d 2 + c 2 d\) 
i ((c 2 di + d\c 2 ) - ( c[d 2 + d 2 c i)) 


V 

(c*di + d 2 c 2 ) — {d\ci + c 2 d 2 ) 

/ 


(5.24) 


(5.25) 


6 Reference frame for Dirac particles 

In this section, we will construct a measurement frame on purely geometric criteria 
and specify the communication protocol. We seek a construction that is coordinate 
independent and as closely tied to the symmetries of Kerr geometry as possible. In 
the photonic case, one could use the principal null directions of the Weyl tensor 
to define a pair of basis vectors for the plane of polarization. Unfortunately, that 
construction cannot serve us in the case of massive spin-1 particles since the spin 
vector cannot in general be confined to a 2-plane, as was the case with the polarization 
vector of a photon. 

Recall that U denotes the 4-velocity of the Dirac particle. We first define a 
volume form for the 3-space (U) 1 - by 

n-=*u\ (6.1) 

where * denotes the Hodge duality operator. Since 

U b := (AE)~5 _ Rhu 1 _ A^Da; 2 - A^a; 3 ) , (6.2) 

the volume form Q,, defined by (16.11) . is given 

(AS) 5 f2 = -A^Q^w 0 A w 1 A a; 2 + A^Bw 0 A w 1 A w 3 

—Rico 0 A a> 2 A a; 3 + Pa; 1 A a; 2 A a; 3 . (6.3) 
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The volume form Q, given by (16.311 should be a monomial. This is indeed the case if 
we consider 1-forms restricted to (C/)- 1 . Let a ! := for i = i, 2, 3, where is 
the parallel propagated frame constructed in Section [2j Then, it follows from (16.3[i 
that 

SI = a* A ct 2 A a' 5 . (6.4) 

Consider the symmetric frame dual to the coframe defined by ( 12 . 1011 . Since one 
of the basis 4-vectors, (e 0 ), is timelike and three (ej,i = 1,2,3) are spacelike, there 
is still an ambiguity in labeling the spacelike basis 4-vectors corresponding to six 
permutations. We show how Carter observers can use geometric criteria to eliminate 
this ambiguity and agree on the labeling of the symmetric frame. 

The symmetric linear combination of the principal null directions of the Weyl 
tensor, 2 “ 2 (£ + n), is timelike, while their difference 2“ 2 (£ — n ) is spacelike. This is 
a geometrically privileged vector that we can label e\. Since the observers’ 4-velocity 
vector e 0 := 2 “z (£ + n ) is known, they can identify their spacelike 4-acceleration 
vector, whose components in the symmetric frame are given by 

V eo e 0 = ( 0 , A 2 ((r 2 — a 2 cos 2 fl) M — ra 2 sin 2 •$) ,0, — a 2 cos sin j . (6.5) 

That is, of the two remaining vectors of the symmetric frame to be labeled, the 
acceleration vector, given by (16.5ji . is orthogonal to only one of them, which we label 
e 2 - The last remaining basis 4-vector is labeled e 3 . Thus, there is no ambiguity in 
labeling the indices of the symmetric frame. 

We now obtain a basis for ([/) _L by contracting the spacelike basis 4-vectors with 
n in an increasing sequence (i < j ) and performing a Gram-Schmidt orthonormal¬ 
ization. One obtains 

W 

Y 

Z 

where 

g := P 2 — A(k — a 2 cos 2 $), p := P 2 — AO, (6.9) 

Definition 6.1 (Reference 3-frame). The set of three spacelike 4-vectors, {X, Y , Z}, 
given by (I 6 . 6 [) ~ (j 6 . 8 p . constitutes an orthonormal basis for ( U that we call the 
reference 3-frame. 


:= p-5 f_A50§,O,O,—P) , (6.6) 

:= (A^PD, 0,p, A©0^ , (6.7) 

:= (ASp)“2 (-FR^,-g, —DA^i?T—A 5 O 2 R 2 j : (6.8) 
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We are now in a position to specify our communication protocol. Let Alice and 
Bob be two Carter observers located in Kerr geometry. In order to communicate 
with Bob, Alice polarizes a massive spin-| at event xa by choosing a unit-norm 3 

component vector W in the reference 3-frame { X , Y, Z } at event xa, and launches 
it on a timelike geodesic y(r) that intersects with Bob’s worldline. Bob intercepts 
the particle and measures its spin vector by projecting it onto the reference 3-frame 
{X, Y, Z} at event x B . The timelike geodesic 7 must satisfy 

7(h)) = x A , and 7(71) = x B , (6.10) 

for some 77 > r 0 . We may without loss of generality set r 0 = 0, and supress the 
subscript for 77 . Thus, 7(0) = xa and 7(r) = x B . 

Armed with the definitions of the reference 3-frame and the communication pro¬ 
tocol, we can define the geometrically-induced precession of a massive spin-1 particle 
as follows. 

Definition 6.2 (Geometrically-induced precession of the spin vector). The preces¬ 
sion of the spin of the Dirac particle is given by the proper time dependent rotation 
A(r) G 50(3) of the reference 3-frame {X,Y,Zj with respect to the 3-frame Lp) 
that is parallel propagated along U and also spans (t/) -1 . More precisely, we have 

At(r) := L ! (t) (r)L (t >,(0), (6.11) 

where L^ k K(0) is the change-of-basis matrix from the reference 3-frame to the parallel 

propagated 3-frame L( fc )(0) at event xa, whereas L*^(r) is the change-of-basis matrix 
from the parallel propagated 3-frame L^)(t) to the reference 3-frame at event x B . 

Remark 6.3. In equation (16. lip , the hatted indices refer to the reference 3-frame 
and the bracketed indices refer to the parallel propagated 3-frame. We are thus 
identifying the tangent spaces along 7 using the parallel propagated frame. Under 
this identification of tangent spaces, the matrix Ah(r), defined by (16.111) . rotates 
Alice’s reference 3-frame to Bob’s reference 3-frame. Thus, 

W B = A (r)W A - (6.12) 

In the next section, we obtain an exact expression for the rotation A(r), and 
obtain a curvature invariant for the curve traced out by the spin vector on the unit 
sphere under the action of A(r). 
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7 Geometrically-induced rotation of the spin vec¬ 
tor 


Direct computation yields 


A(r) := 


( rDP+g cos i?V _RQ 
x/ftXp 

rg& 2 —acostTOPit^ 
VkS gp 

a cos i?P 


l l 

_k ^ sin \ SP© 2 +ci cos \ 

a? 2 (re+r 2 )v / KE/9 

ft 2 sm\EBP 2 +C 2 cos x 

(/c+r 2 )-\//cUgp 

I cosyrP—sinyi?2 

W 2 --- ■= --- 

y/M> 


1 1 
cos y£PQ2 —ci siny 

ZD 2 (re+r 2 ) v'kSp 

—ft 2 cos \EBP 2 +C 2 sin y 

1x7^ (/■c-(-7’ 2 ) v / k;U£P 

A sinyrP-l-/c2 cosy.R2 

tU2 -^-i—T=--- 

V*Q 


(7.1) 


where 


Ci := (« + r 2 ) yrwVRQ — acos^OPj , C 2 := (k + r 2 ) ^acos'df)© 2 + rtuDPi? 2 

(7.2) 

and 07 is defined by equation (12.331) . while g and p were defined in equation (16.9I) . 
In the matrix on the RHS of equation (17. lj) . all entries depends on proper time r 
through r(r) and $(r). We verify that A T A = / and det A = 1. Thus, A e 50(3). 

Let Wo be the initial position of the spin vector prepared by Alice. At proper 
time r > 0, the position vector is given by 


W(t) = A(r)W 0 , 


(7.3) 


where A(r) is given by (17.11) . 

The differential equation satisfied by the Darboux frame (^W, W , W x W ^ for 
a curve on the unit sphere can be written as \G^n\ 


l w \ 
w 

\W xW J 


0 1 

° \ 

( w 
/ 

\ 

-1 0 

K 

w 


c 

•-Sc 

1 

G 

o / 

->■ ->• 

) 

y 

\W xW 


(7.4) 


where the prime denotes derivatives taken with respect to the arclength parameter 
s and k g (s) is the spherical curvature Jf| Any two curves confined to the unit sphere 

3 The spherical curvature k g is usually called geodesic curvature. But we shall not that termi¬ 
nology to avoid confusion between geodesics on the sphere and geodesics of Kerr geometry. 
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with the same spherical curvature are related to each other by a constant rotation 
and converselyo 

From equation 07.41) . we obtain the following straightforward formula for the 
spherical curvature of curves on the unit sphere parametrized by the arclength s. 


k g (s) 


W ■ W xW 


w ■ w 


(7.5) 


Note that our position vector (17.31) is parametrized by proper time, r, not by ar¬ 
clength s of the curve traced out by the position vector W. We therefore amend 
equation (17.5ft to obtain the following formula for the spherical curvature invariant 
in terms of proper time r. 


M r ) 




W x W) 



(7.6) 


where || • || denotes the Euclidean norm. The spherical curvature k g defined by the 
differential equation (17.4j) and given by formula (I7.6[) is independent of the choice of 

initial position vector Wo- Another choice of initial position vector merely induces 
a constant rotation of the curve that leaves the spherical curvature k g invariant. We 
are thus free to choose the initial position vector so as to minimize computational 
complexity. We fix the initial position vector at r = 0 to be given in the reference 
frame (j 6 . 6 [) - (l 6 . 8 [) by 

' 1 


W n = 


(7.7) 


At proper time r > 0, the spin vector is then given simply by 


/ 


W{t) = A(r)W 0 = 


rDP+q cos 
_ VkE p 

rg\/0~a cos i?DP\7R 
\/kE gp 
a cos i?P 


\ 


y/KQ 


(7.8) 


We are now in a position to prove the following result. 

4 More precisely, if (a;(s)) se ( 0 b) an d (/^( s ))se(o 6) are ^ w0 curves on the un it sphere parametrized 
by arc length s, then there exists A G SO( 3) such that a(s ) = A/3(s),Vs G (a, 6) if and only if 
kg(s) = kP(s),Vs G (a, b). 
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Proposition 7.1. There is no spin precession for Dirac particles confined to the 
equatorial plane, Eq := {—oo < t < +oo,r + < r < +oo,d = 7t/2, 0 < ip < 2i r} ; of 
Kerr geometry. 


Proof. For timelike geodesics confined to the equatorial plane, we have $ — |, g — 
P 2 — k A, w = pfipp and k = (aE — <f>) 2 . Thus, 


W(t) = 



(7.9) 


for arbitrary proper time r. 


□ 


Corollary 7.2. There is no spin precession for Dirac particles in Schwarzs child 
geometry. 

Proof. Since Schwarzschild geometry is spherically symmetric, geodesics are confined 
to planes through the origin |Cha92| . Any plane through the origin in Schwarzschild 
geometry is isometric to the equatorial plane of a degenerate Kerr solution with 
a = 0. □ 


In Section [8j we will discuss these results as well as present some plots for quali¬ 
tative analysis. 


8 Summary and discussion 

In order to isolate the geometrically-induced precession of the spin of Dirac particles 
we started with the Dirac equation. In general, there is no way to associate a solution 
of the Dirac equation (13.111 to a classical trajectory; that is, a timelike geodesic along 
which the nonzero rest mass, spin-I test particle propagates. Following the insight 
of |Pau32j . we deploy the semiclassical ansatz for the Dirac spinor. The classical 
limit {h —y 0) for the Dirac equation with the semiclassical ansatz is equivalent to 
the Hamilton-Jacobi equation for spinless particles. This yields the desired geodesic 
that we postulate is the classical trajectory of the Dirac particle. 

In order to define the spin vector, we appealed to the decomposition of the Dirac 
current due to Gordon (Theorem 13.IF The spin vector defined in terms of the 
conserved polarization current in the Gordon decomposition (Definition 13. 2\i . corre¬ 
sponds to the Foldy-Wouthuysen mean-spin operator |KRH01j . Following the dis¬ 
cussion in Section |3] we can be confident that this definition of the spin vector, first 
proposed in |Aud81j . is a good one. 
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Extending the result of [AudSlj . we showed that the spin vector defined in terms 
of the conserved polarization current in the Gordon decomposition, is parallel prop¬ 
agated along the geodesic obtained from the semiclassical ansatz to 0(h 2 ) in Kerr 
geometry (Theorem [5]). We were somewhat surprised by the result as we were ex¬ 
pecting a more involved propagation law for Dirac particles. 

We showed how observers located in the outer geometry of Kerr black holes may 
set up reference 3-frames in terms of locally-available, purely geometric information 
(Definition 16.11) . We were thus able to define the geometrically-induced precession 
of the spin vector of Dirac particles propagating in the outer geometry of Kerr black 
holes (Definition 16. 2 j) . Our geometrically-motivated strategy allowed us to obtain a 
compact expression for the geometrically-induced spin precession for Dirac particles. 
We showed how the geometrically-induced precession of the spin of Dirac particles 
is determined by the rotation of the parallel-propagated frame with respect to the 
reference frame that the observers use to measure the spin of the test particles. We 
further showed that the geometrically-induced spin precession of a Dirac particle can 
be invariantly represented by the spherical curvature of the curve traced out by the 
spin vector on the unit sphere. 

We showed how there is no spin precession for Dirac particles confined to the 
equatorial plane of Kerr geometry (Proposition 17.11) . The significance of this result 
comes from the fact that many authors restrict attention to the equatorial plane in 
order to simplify computations (e.g., [RPGClll ISA10] ). Note that the non-triviality 
of the results obtained in [RPGGlll ISAIO] is the result of the non-zero acceleration 
of their chosen test particles. 

Finally, we showed how, as an immediate consequence of Proposition 17.11 spin 
precession is trivial for Dirac particles in Schwarschild geometry (Corollary 17.21) . 

In order to qualitatively analyze the expression we have obtained for the spin pre¬ 
cession of Dirac particles, we present some plots. In each of 3 sets of figures, the first 
figure (a) shows the orbital behaviour of the timelike geodesic with (r(s), <p(s)) as 
polar coordinates, the second figure (b) depicts the same orbits in three dimensions 
with spherical coordinates (r(s), i9(s), ip(s)), and the last figure (c) depicts the spher¬ 
ical curvature k g as a function of propertime r. Figures 1 and 3 present co-rotating 
orbits, while Figure 2 presents a counter-rotating orbit. 
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Figure 1: A counter-rotating orbit with E = 2, = 3, k = 12 and initial data 

r(0) = 20,79(0) = 1.57, and y>(0) = 0. 



(a) The orbit in polar coordinates (x = rcos(p, y = r sin tp). 



(b) The orbit in 3D spherical coordinates {x = r cos ip sin i9, y = r sin (p sin -d, * = rcosd). 



(c) The spherical curvature k g as a function of propertime. 










Figure 2: A counter-rotating orbit with E = 1.004, = —4, k = 60 and initial data 

r(0) = 20,79(0) = 1.57, and y>(0) = 0. 



(a) The orbit in polar coordinates (x = rcos(p, y = r sin tp). 



(b) The orbit in 3D spherical coordinates {x = r cos ip sin i9, y = r sin (p sin -d, * = rcosd). 



(c) The spherical curvature k g as a function of propertime. 
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Figure 3: A counter-rotating orbit with E = 1.004, $ = 4, k = 16 and initial data 
r(0) = 20,0(0) = 1.57, and y>(0) = 0. 



(a) The orbit in polar coordinates (x = rcos(p, y = r sin tp). 



(b) The orbit in 3D spherical coordinates {x = r cos ip sin i9, y = r sin (p sin 0, z = r cos??). 



(c) The spherical curvature k g as a function of propertime. 
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